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Applying a Keldysh Green's function method it is 
shown that hot electrons injected from a STM-tip into a 
CoSi2/Si(lll) system form a highly focused beam due to 
the silicide band structure. This explains the atomic res- 
olution obtained in recent Ballistic Electron Emission Mi- 
croscopy (BEEM) experiments. Localized surface states in 
the (2 x l)-reconstruction are found to be responsible for the 
also reported anticorrugation of the BEEM current. These re- 
sults clearly demonstrate the importance of bulk and surface 
band structure effects for a detailed understanding of BEEM 
data. 

PACS numbers: 61.16.Ch, 72.10.Bg, 73. 20. At 

Ballistic Electron Emission Microscopy (BEEM), and 
its spectroscopic counterpart (BEES), are powerful tech- 
niques invented for exploring the electronic properties of 
metal-semiconductor (M-S) interfaces Thin metal- 
lic films are deposited on different semiconductor mate- 
rials and the BEEM current, i.e. the current arriving 
at the semiconductor after injection into the metal sur- 
face from a Scanning Tunneling Microscope (STM) tip, 
is measured as a function of the tip- metal voltage 
The interpretation of these experiments is based on a 
three-step model: (i) first, electrons are injected from the 
tip into the metal (tunneling); (ii) then, electrons prop- 
agate through the film suffering collisions with different 
quasiparticles (transport), and (iii) finally, electrons over- 
come the Schottky barrier and enter into the semiconduc- 
tor (matching of metal and semiconductor wavefunctions 
across the interface). The difficulty in analyzing exper- 
imental BEEM data stems from the strong influence of 
all three steps requiring a careful theoretical modeling to 
avoid spurious correlations between the parameters in- 
volved. Recently, it has been shown that the electronic 
band structure of the metal, which had been completely 
neglected in earlier free electron models, plays a crucial 
role in this regards || . 

Recent experimental BEEM investigations on metal- 
lic silicide films deposited on Si show (i) an atomic scale 
resolution of the M-S interface and (ii) a striking de- 
pendence of the interface BEEM current on the silicide 
surface topography ||. Dislocations and point defects 
at the interface were well visible giving direct access to 
its quality. This is rather important in view of the in- 



terface's role in building the Schottky barrier or with re- 
spect to the growth mode of silicides which are promising 
materials for microelectronic applications pj. Equally 
important is the quantitative understanding of how the 
obtained atomic scale resolution is produced and why the 
BEEM current is related to the surface topography. In 
this paper, we show that the high lateral resolution is 
caused by the silicide's band structure, which in the case 
of CoSi2/Si(lll), on which we concentrate, makes the 
electrons focus in the <111> direction. This tells that the 
experimentally observed focusing is an intrinsic feature of 
such films that might be exploited in future applications. 
Additionally, the introduction of the appropriate surface 
electronic structure explains the BEEM current depen- 
dence on the tip position mainly as a result of the weight 
of localized surface states on the reconstructed surface. 

We use a full quantum-mechanical description of the 
BEEM problem based on a Keldysh Green's function 
method 0. This formalism presents the important ad- 
vantage over standard E-space Monte-Carlo approaches 
of yielding an appropriate description of the electronic 
band structure. Moreover, inelastic effects associated 
with electron-electron interactions are also included in 
our method by adding a positive imaginary part to the 
energy of the electron. In order to analyze the first two 
steps of the BEEM process, we choose a local orbital 
basis for the description of the electronic structure of 
the tip and sample and the coupling between them. In 
particular, for CoSi2 we use a slight modification of the 
tight-binding parameters given in Q , that accurately re- 
produce the band structure of this silicide around the 
Fermi level Q. For the analysis of the interaction be- 
tween tip and sample we assume that only the last atom 
in the tip (0) is connected to the sample. Hence, we ex- 
press the coupling tip-sample in terms of a set of hopping 
matrices Xb m , that link the tip atom (0) with the atom 
(m) in the sample surface. For each To m , a WKB de- 
rived exponential damping is applied, valid because the 
tip-sample distance in BEEM is rather large. 

Being interested in understanding the observed nano- 
metric spatial resolution of this technique, we first an- 
alyze currents in real space. Within our formalism, the 
current between two sites i and j in the metal can be 
obtained from the following formula 
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J/(S, ky) = ^3Tr 



where is the hopping matrix linking local orbitals of 
both sites (i and j), and the trace denotes summation 
over these orbitals. g^ m (E) is the retarded Green's func- 
tion for the surface decoupled from the tip. This function 
describes the propagation of an electron between atoms 
j and m inside the metal, including the effect of the sur- 
face. Atom to in Eq. (1) is coupled to the tip atom by 
a hopping matrix T m Q, and poo(E) is the density of states 
matrix at the tip atom. The advanced Green's function, 
g£i (E) , describes the electron propagation from an atom 
7i at the surface down to the atom i, closing the loop to 
give the current between atoms j and i. The summa- 
tion runs over all tunneling active atoms in the sample 
surface, m and n. The energy integration is performed 
between the Schottky barrier (eV , assumed to be 0.66 
eV j|) and the applied voltage (eV). However, due to 
the exponential energy dependence of the coupling ma- 
trices, T m o,To n , the integrand is a strongly increasing 
function with energy, so that already the contribution at 
the highest energy (eV) provides the dominant fraction 
of the elastic BEEM current in the near threshold region. 
To elucidate the physics behind the observed effects, the 
presented results will therefore be restricted to this high- 
est energy. 
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FIG. 1. Current distribution in a Si2 layer parallel to the 
surface after propagation through 30A CoSi2(lll) film. In- 
jection from the tip at 1.5eV occured in the center of the 
shown plane (white X), where the maximum current propa- 
gating in a focused beam along the <111> direction can still 
be found. The linear gray scale indicates current intensity at 
each atomic site: black maximum to white zero current. The 
right hand panel displays a cut through the focused beam in 
<211> direction from which a FWHM of 8.9A can be derived. 

With Eq. (0), the elastic propagation of electrons in 
real space from the tip down to the M-S interface can 
be followed. In order to obtain the final BEEM current, 
we further need to calculate the momentum distribution 
of the electrons that reach the M-S interface, Ji(E, kn). 
This momentum distribution can be expressed as B : 
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where in this case g^(E, k||) is the retarded Green's func- 
tion for the unperturbed metal, linking layer c (the metal 
layer at the M-S interface) and the surface layer 1 which 
is connected to the tip by a hopping matrix Tio(kii). 
Tf, c (k||) is the hopping matrix connecting all upper lay- 
ers b with the interface layer c and finally ^(.E, k||) is 
the advanced Green's function linking the surface layer 
with layer b. These advanced and retarded Green's func- 
tions and the ones appearing in Eq. ([!]) can be readily 
computed using renormalization group techniques Efl] , 




FIG. 2. Electronic current distribution in the 2D interface 
Brillouin zone, Ji(E,k\\), evaluated at 1.5eV after 30A film 
propagation. The current intensity is drawn with a linear 
gray scale, black representing maximum current. Also shown 
are the ellipsoids defining available states in the semiconduc- 
tor below 1.5 eV. The inset contains the constant energy sur- 
face sheet mainly responsible for the current propagation: the 
shaded flat terraces point in <llt> direction and correspond 
to the dark areas of the 2D current distribution. 

First of all, we apply our formalism to the case of a 
CoSi2(lll)-(lxl)Si-rich surface terminated film [jllj and 
analyze the propagation of electrons from the tip to the 
CoSi2/Si(lll) interface in real space. In the <111> di- 
rection, this metal may be characterized by a stacking 
sequence of Sil-Co-Si2 trilayers (cf. Fig. 3a) with the 
interface to Si mainly formed below a Si2-type layer [pd] . 
Fig. ^ shows the current distribution on every atom in 
such a Si2 layer 30A below the surface as to compare 
with experiments performed on films of equal width [0. 
The prominent effect we deduce from this figure is that 
the electrons injected into the silicide are focused inside a 
very narrow beam propagating perpendicular to the film. 
The right hand panel in Fig. 1 shows the intensity in real 
space along a line in <211> direction through the cen- 
ter of the beam: the obtained FWHM of 8. 9 A compares 
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very well with the resolution of « 10. OA with which in- 
terface point defects could be resolved experimentally in 
such films B. This, up to now, highest achieved spatial 
resolution with BEEM had been impossible to explain 
assuming free electron propagation inside the metal, pre- 
dicting beam widths of 25 A for the same distance. 

The electron focalization is due to the particular shape 
of the constant energy surface sheet responsible for the 
major current propagation (see inset of Fig. |J); it can be 
shown that between Ep and Ep + 2.5 eV these sheets are 
practically the same except by a uniform shrinkage that 
increases linearly with energy. The shaded regions are 
nearly flat terraces perpendicular to the <L11> direction, 
and act as a kind of "condenser lens" on the electron 
beam, keeping the electrons with corresponding k-vector 
propagating along the <111> direction ||. This reason- 
ing is complementary to the current distribution we have 
calculated in k|| -space using Eq. 2 and shown in Fig. [|. 
The three dark regions of the 2D Brillouin zone where 
the k||-currcnt is mainly concentrated correspond to the 
flat areas of the constant energy surface. 

The onset of BEES I(V) characteristics is linked to 
the Schottky barrier height between the metal and the 
semiconductor. In our calculations that onset appears at 
0.9eV, 0.24eV larger than the Schottky barrier height 
commonly accepted for the CoSi2-Si interface fll2| , |l3|] . 
This is related to the assumed k|| -conservation and to 
the absence of states in the metal matching the conduc- 
tion band minima in the semiconductor. The same de- 
layed onset has been obtained by Stiles and Hamann [jl4| , 
and it has been argued 13 13] that a smaller onset can 



appear if a non ky-conserving scattering process is oper- 
ative for the injected electrons at the silicide-silicon in- 
terface. Indeed, the results reported by these authors for 
CoSi2/Si(lll) seem to point out that the effect of such 
processes is to modify only slightly the BEEM-current 
beyond 0.9 eV, but is enough to yield the appropriate 
M-S barrier height at w 0.66 eV. Therefore, for energies 
larger than 0.9 eV, current injection conserving k|| dom- 
inates the spectra, as expected intuitively from the good 
matching between the Si and CoSi2 lattices, and our the- 
ory applies. We should also mention, however, that non 
k|| -conserving processes must play an important role for 
the BEEM contrast of defect images at the M-S interface 
due to the nanometric size of the electron beam. 

In our next step we consider the case of the (2 x 1) 
surface structure of CoSi2/Si(lll). Stalder et al. jl6| re- 
ported this Si-rich reconstruction for CoSi2(lll) with a 
geometry very similar to Pandey's 7r-bonded chain model 
p7| . Fig. 3a shows a sideview of this surface geometry 
with its topmost Si-bilayer reconstructed in alternating 
high and low chains. We have analyzed how the geome- 
try of the reconstructed (2x1) surface modifies the elec- 
tron focalization discussed above for the (lxl) surface. 
The main effect of the reconstruction is to broaden the 
FWHM of the focused beam to 13.6A. This effect, that we 



associate to a larger area of the surface unit cell where the 
tunneling electrons are injected, has also been observed 
experimentally by Sirringhaus et al . 

A very interesting result observed for this reconstruc- 
tion is that in the constant-current STM mode, the 
BEEM image of the interface reflects the atomic surface 
periodicity, but out of phase with the topographic cor- 
rugation H . This BEEM anticorrugation has previously 
been attributed to atomic-scale variations of the energy 
tunneling distribution of the injected electrons Q. In or- 
der to analyze these results, we calculate the current that 
reaches the M-S interface as a function of the tip posi- 
tion for a constant tunneling current (1 nA). In general, 
to compute the BEEM current injected into the semicon- 
ductor we would need to use a transmission coefficient, 
T(E, k||), determined from the matching of states at the 
interface. However, this is not necessary to study the 
particular dependence of the BEEM current on the tip 
position, as T is independent of the tunneling injection. 
Moreover, we have found that the k|| distribution of the 
current is nearly the same for all the different positions 
of the tip, in accordance with the conclusions raised in 
Therefore, to study the effect that the surface recon- 
struction introduces in the BEEM current we can simply 
analyze the total current reaching the interface. Fig. 3b 
shows that this quantity presents anticorrugation with 
respect to the one found in the surface reconstruction. 
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FIG. 3. a) Sideview of the (2 x 1) surface reconstruction 
with tip positions 1 (high chain injection) and 2 (low chain 
injection), b) M-S interface current at 1.5eV after 30A film 
propagation as a function of the tip position (higher curve, 
left scale in pA). The lower curve gives a schematic surface 
topography along the scan line (right scale in A). 

To understand the physics behind this effect, we have 
studied the injected current along the different metal lay- 
ers. Fig. 4 shows our results for the tip located either on 
the highest or on the lowest position on the reconstructed 
surface (points 1 and 2 of Fig. 3a). Two important con- 
clusions can be drawn from our results: first, the injected 
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current along the metal layers is damped by the introduc- 
tion of an imaginary component for the energy, E + in, 
that simulates the electron-electron scattering processes 
(in our case we have used r\ = 0.05eV, that yields at- 
tenuation lengths in accordance with experimental data 
pi!); only at long distances this damping results in an 
exponential behaviour for the current. Second, at short 
distances the current presents a faster decrease associ- 
ated with the injection of electrons along surface states 
channels. As Fig. 4 shows, in the (2 x 1) reconstruction 
the current decreases by 65% after the electrons cross 
the first two Si-layers and the first CoSi2 trilayer (where 
surface states are mainly localized). This is the effect 
of having the injected electrons propagating also along 
the surface bands, departing in this way from the bulk 
states channels that contribute to the current propagat- 
ing across the metal layers. 
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FIG. 4. Current across the first 30A of the silicide film at 
1.5eV after injection on a high chain (solid line) or on a low 
chain (dotted line). The diamonds at the bottom indicate the 
trilayer Sil-Co-Si2 sequence of CoSi2 in the (111) direction, 
with the final reconstructed Si chain-bilayer. The inset shows 
the surface density of states projected on the high and low 
chain atoms in the energy region important for BEEM. 

The anticorrugation obtained in Fig. 3b for the total 
current arriving at the M-S interface can be understood 
in terms of this role played by the surface states. These 
states have larger weights on the atoms of the higher 
chains than on the low chain ones as we can see in the 
surface density of states shown in the inset of Fig. 4. 
Therefore, for those tip positions in which electrons are 
injected predominantly into high chain atoms, there is 
a larger probability for having those electrons propagat- 
ing along surface states channels than the one obtained 
for injection on the low chain atoms. Consequently, in 
the usual experimental constant-current STM mode, less 
current crossing the metal layers and reaching the inter- 
face remains. This difference with respect to the current 
injection, which stays almost constant after crossing the 



first CoSi2 trilayer (5A), explains why the BEEM image 
shows anticorrugation with respect to the surface topog- 
raphy. Note, however, that the absolute order of mag- 
nitude of this effect depends strongly with energy: as 
shown in Fig. 3b, the anticorrugation contrast is 25% for 
1.5 eV and a lower contrast is obtained for larger volt- 
ages. This dependence is related to the fact that surface 
states are concentrated rather close to the Fermi level. 

In conclusion, we have presented a theoretical analy- 
sis of the propagation of an electron beam injected in a 
CoSi2(lll) crystal using a STM tip. Our results show 
conclusively that the silicide electronic band structure 
plays a central role in the focalization of the electron 
beam. This behaviour and the specific k|| -contribution to 
the current have been associated with flat terraces of the 
constant energy surface producing a condenser lens effect 
on the electron propagation. Our results explain the high 
resolution observed in real space for BEEM experiments 
performed on CoSi2/Si(lll) interfaces. Additionally, we 
have also shown how the BEEM current can map out 
the silicide surface reconstruction due to the role played 
by the localized CoSi2(lll)-(2xl) surface states on the 
current injected from the tip. 
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